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Abstract 

Structural, electronic, and optical properties for the cubic, tetragonal, and monoclinic crystalline 
phases of Zr02, as derived from ab initio full-relativistic calculations, arc presented. The electronic 
structure calculations were carried out by means of the all-electron full potential linear augmented 
plane wave method, within the framework of the density functional theory and the local density 
approximation. The calculated carrier effective masses are shown to be highly anisotropic. The 
results obtained for the real and imaginary parts of the dielectric function, the reflectivity, and 
the refraction index, show good agreement with the available experimental results. In order to 
obtain the static dielectric constant of Zr02, wc added to the electronic part, the optical phonons 
contribution, which leads to values of ei(0) ~ 29.5,26.2,21.9, respectively along the xx,yy, and zz 
directions, for the monoclinic phase, in excellent accordance with experiment. Relativistic effects, 
including the spin-orbit interaction, are demonstrated to be important for a better evaluation of the 
effective mass values, and in the detailed structure of the frequency dependent complex dielectric 
function. 
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I. INTRODUCTION 



Zirconia (Zr02) is a material of great technological potential importance due to its out- 
standing mechanical and electrical properties, high dielectric constant and wide band gap. 
Among the Zr02 applications, are: its use as gas sensors, solid fuel cells, high durability 
coating, catalytic agents, etc. In recent years, its gap {Eg ~ 6 eV) and dielectric properties 
{e ~ 25) suggested its potential to replace Si02 in advanced metal oxide semiconductor 
devices (MOS) in gate stack, dynamic access memory devices, and optical applications [l-5|. 
Moreover, it has a large band-offset in direct contact with Si and good thermal stability. 
These attracting properties has led zirconia based oxides to be widely studied in recent years 

Zr02 presents polymorphism with monoclinic, tetragonal, and cubic phases. The mono- 
clinic phase (m-Zr02) of undoped Zr02 is thermodynamically stable at temperatures below 
1170 °C. The tetragonal phase structure (t-Zr02) arises when heating Zr02, being stable 
between 1170 °C and 2370 °C, as well as above this; through its melting point (2706 °C), 
the cubic phase (c-Zr02) is then observed. The temperature in which the tetragonal to cubic 
transformation occurs can be lowered, by the addition of solutes such as MgO, CaO, Y2O3, 
allowing the achievement of the stabilized c-phase even at room temperature. 

Previous investigations of Zr02, based on first-principles electronic structure calculations, 
have been reported in the literature j9-16|. Most of them dealt with the structural and elec- 
tronic properties of Zr02, in its cubic, tetragonal, and monoclinic phases IJQ, Q. The 
zirconia optical properties, or more specifically its complex dielectric function has been ad- 
dressed by French et al (electronic frequency regime only) by carrying out ab initio electronic 
structure calculations via the orthogonalized linear-combination-of-atomic-orbitals method 
3. and ,„o..e .cent,, bv Zhao and Va„de..bUt Q and by Rignanese et by mean, of 



ah initio pseudopotential calculations. However, in the work of Ref. [17[ only the Zr02 cubic 
and tetragonal structural phases have been considered, while in that of Ref. 1^, although 
in the work all three phases (cubic, tetragonal, and monoclinic) are studied, only the lattice 
contributions to the dielectric tensor have been addressed so far. Moreover, in these studies, 
relativistic effects including the spin-orbit interaction were not taken into account in the 
calculations, neither have been provided values for the Zr02 carrier effective masses, which 
are fundamental to model the tunnelling currents through floating gate memory devices that 
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use, e.g., Zr02 as the tunnelling oxide. 

In this work, we present results of state-of-the-art first-principles calculations for the 
structural, electronic, and optical properties of Zr02 in its monoclinic, tetragonal, and cubic 
phases. From the band structures, the Zr02 conduction- and valence-band effective masses 
have been obtained, in the most important symmetry axis of the Brillouin zone. The fuU- 
relativistic effects, properly taken into account in the calculations, are demonstrated to be 
relevant for a more precise evaluation of its electron- and hole-effective-mass values, as well as 
in the detailed structure of the electronic contribution to its frequency dependent dielectric 
function, As a consequence of the role played by the inclusion of full-relativistic 

effects in the calculations, the obtained changes in the reflectivity spectrum and refraction 
index are also significant. 

In order to realistically obtain the value for the static dielectric constant ei{0), we also cal- 
culated the lattice contributions, mostly arising from the transverse optical phonon modes, 
to the low-frequency dielectric function, for the different phases. The resulting calculated 
static dielectric constant for m-Zr02 is in excellent agreement with the values reported in 
the literature. Moreover, the results obtained for the reflectivity and refraction index for its 
monoclinic phase are in good accordance with the available experimental data. 



II. CALCULATION DETAILS 



The electronic structure calculations were carried out using the density-functional the- 
ory in the local-density approach (DFT-LDA) by means of the ab initio full-potential self- 
consistent linear augmented plane wave (FLAPW) method (Wien2k code 18|). The gen- 
eralized gradient approximation (GGA) was adopted for the exchange-correlation potential 
jiol l. A similar approach has been recently applied to study the structural, electronic, and 
optical properties of cubic SrTiOa and Hf02 2Q|, |2l|. The Zirconium 4s-,4p-,4d-,and 5s, and 
the Oxigen 2s, 2p electrons were treated as part of the valence states. The cutoff angular 
momentum was / = 10 for wave functions, and / = 5 for charge densities and potentials 
inside the spheres. Muffin-tin sphere radii were assumed equal for both Zr and O. The 
number of symmetrized k-points, used as input for the self-consistent charge density de- 
termination, was 47, 75, and 100-k points in the irreducible symmetry wedge of the BZ, 
respectively for the cubic, tetragonal, and monoclinic phases. The value RKMAX = 9 was 
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adopted for the three different phases. With these assumptions, the energy bands were con- 
verged within 10~^ eV, and the total energy within 10~^ eV. The core electron states were 
treated full-relativistically, whereas the valence states were treated both, non-relativistic 
and full-relativistically, i.e. including the spin-orbit interaction, besides the Darwin and 
mass-velocity corrections. 

In addition, we have also used the DFT-LDA, with the plane-wave description of the 



wavefunctions, and the pseudopotential method {Abinit code 22|) for the evaluation of the 
phonon contributions to the real part of the dielectric function {e^i^^{w)), as will be discussed 
later. With this approach, we are able to unveil and estimate the total value of the static 
dielectric constant of Zr02 in the form £1(0) = ef^'^'^{0) + £:5"**(0), for the three phases cubic, 
tetragonal, and monoclinic. 

III. RESULTS AND DISCUSSIONS 
A. Structural properties 

The cubic Zr02 phase (c-Zr02) belongs to the space group Fm3m with crystalline struc- 
ture of fluorite, i.e. face-centered cubic lattice with three atoms at the basis, being one 
Zr atom at the position (0,0,0) and the oxygen atoms at the position (1/4, 1/4, 1/4), for 
which the BZ is the usual 14-faces polyhedra. We have obtained through a total energy 
minimization process, a lattice constant value of the cubic phase a = 5.139 A, slightly larger 
than that calculated by Kralik et al [l^, using the pseudopotential method and the local 
density approximation (LDA), and closer to the experimental value of 5.09 A 13|. The value 
obtained by us confirms that the GGA approximation leads to lattice constant values which 
are slightly overestimated in relation to the experimental value. 

The tetragonal phase (t-Zr02) is described by the PA2/nmc group. As shown by Zhao 



and Vanderbilt the t-Zr02 phase is obtained from the cubic one, by properly performing 
alternate displacements Ac of the oxygen atoms in the cubic structure. Our optimization 
process consisted in minimizing the total energy and forces, until "converged" Ac, c, and a 
values have been reached. The obtained values were a = 5.10 A, c = 5.23 A, and the internal 
parameter d^ = Ac/c = 0.05, close to the experimental ones [13j (a = 5.05 A, c = 5.18 A, 



and d^ = 0.0574), as well as close to the theoretical values reported by Kralik at al. [12| 
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(a = 5.04 A and c = 5.10 A and d, = 0.0423). 

For the monoclinic phase, which is described by the P21/c group, we have evaluated the 
total energy for several lattice constant values, in order to have a good description of the 
equilibrium properties for the bulk m-Zr02, in which case, for sake of reducing computational 



time, we use the Abinit code 22|. The results were fitted by the Murnaghan equation of 
state, and we have obtained, for the monoclinic structure, a = 5.12 A, b = 5.16 A, c = 5.33 
A, and e = 99.6 de grees. These results are in accordance with the available experimental 



13|(a = 5.15 A, b = 5.21 A, c = 5.31 A, and 6 = 99.23 degrees), and theoretical results 
15j(a = 5.20 A, 6 = 5.25 A, c = 5.41 A, and 6 = 99.60 degrees). We have used the later 
ones in the band structure calculations performed with Wien2k code, as well as to calculate 
effective mass values and optical properties for the monoclinic phase. 

B. Electronic properties 

Table I presents the band gap energies, together with the valence-to-conduction band 
transition symmetry, for all three phases of Zr02. For the cubic phase, the band gap is 
indirect at X — > F, with an energy of 3.09 eV when a non-relativistic calculation is consid- 
ered, and with energy of 3.30 eV, if the calculations are performed by taking into account 
the full relativistic effects, i.e. the scalar-relativistic and spin-orbit effects. The direct band 
gap, at F, is of 3.61 (3.80) eV within a non-relativistic (fuU-relativistic) calculation. Our 
values for the band gaps are very close to those obtained by Kralik et al. [12I, in which a 
pseudopotential LDA method was used. We recall that, the theoretical values for the band 
gap energy are smaller when compared with the experimental one, due to the well-known 
underestimation of conduction band states energies in ab initio calculations which are per- 



formed within DFT. Experimental values, as reported by French et al. 16|, are found in 
the range of 6.1 to 7.08 eV. 

Due to the spin-orbit interaction, we observe a spin-orbit (so) splitting energy of the 
valence band top (f), at the F-point, which for the cubic phase is A^o(F) = 69 meV. There 
are no experimental data for A^o in Zr02 reported so far. 

For the tetragonal phase, the energy transitions at X — v F, Z — > F, and F — > F are 
very near. An indirect gap of 3.80 (4.01) eV, and a direct one of ~ 3.9 (~ 4.1) eV, were 
obtained from non-relativistic (full-relativistic) calculations. A value of A^^(F) = 9 meV 
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has been observed for the spin-orbit sphtting energy of the valence band top, at F, for t- 
Zr02. The band gap energy values, obtained through the fuU-relativistic calculations, are 
well compared with those obtained by Kralik et al. 12|. The experimental values for t-Zr02 
lie in the range of 5.8 to 6.6 eV, according to the results reported by French et al. jl6|. 

For the monoclinic phase, we have obtained an indirect gap of ~ 3.4 eV (~ 3.6) eV in 
a non-relativistic (full-relativistic) calculation, and a direct band gap of 3.50 eV (3.64) eV, 
so much close to the value of the indirect one, that it is difficult to assure whether the gap 
is direct or indirect. As one goes from the cubic to tetragonal, and then to the monoclinic 
structural phase, the symmetry lowering is responsible for the removal of the degeneracy, 
hence to the appearance of so many energy bands in the monoclinic case. 

Figure 1 shows the non-relativistic results obtained for the band structure, along high 
symmetry directions of the BZ, and the total density of states (TDOS) for the three phases 
of Zr02. The main character of the peaks in the TDOS is also emphasized. The results for 
band structures, as obtained from the full-relativistic calculations are presented in Fig. 2. 
In the upper valence part of the band structure, the 0(2p)-related states are predominant, 
while in its lower part, the 0(2s)- derived states are the most relevant, independently of 
inclusion of the relativistic corrections, for all three phases. The Zr(4p) states appear rather 
deep in the valence band (-26 eV) and, due to the spin-orbit coupling are split into a group 
at -26 and a group at -28 eV, corresponding to j=3/2 and j=l/2 manifolds, respectively. 



C. Carrier effective masses 



The valence- and conduction-band effective masses for the cubic, tetragonal, and mon- 
oclinic phases of Zr02, at relevant symmetry points of the BZ, are presented in Table II 
for several directions. Full- and non-relativistic (see the values in parenthesis) results are 
shown. In order to remain within the region of validity of the parabolic approximation, the 
effective masses were calculated by considering the energy curves within a small region of 
radius 0.5% centered on the extremum of interest. 

As we observe, in both the valence and conduction bands, the carrier effective masses are 
very anisotropic, i.e., show a relevant dependence on the k direction. Moreover, remarkable 
differences arise when the relativistic effects are taken into account. In particular, for the 
cubic phase, relativistic effects are seen to be most important in the T — X direction, for 
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holes and electrons, while these effects are more relevant in the T — Z direction, for the case 
of holes, for the tetragonal phase, and in the T — Y direction, for electrons and holes, for 
the monoclinic phase. So far, there are no reported experimental data in the literature of 
carrier effective masses in Zr02. This represents the first report of effective mass values, as 
calculated for all phases of Zr02, of fundamental importance to the modelling of advanced 
devices based on this oxide. 



D. Optical properties 



The electronic contributions to the real {ei{u))) and imaginary {62(00)) parts of the com- 
plex dielectric function, e{u) = ju) + 62(00), can be obtained from the band structure 
directly through the Wien2k code 18|, ei from €2, and vice- versa, using the Kramers-Kronig 



relations 



23| . Once is known the imaginary part 62(00), from the FLAPW electronic struc- 



ture calculations, the real part ei(oo) is then given by [2 



Other optical properties, such as the reflectivity R(w) and the complex refraction index 
n(ci;) can also be obtained. We took into account that energy gaps are under evaluated by 
DFT [1^, by adopting the following procedure: the dielectric function, and all the related 
optical properties, were obtained by performing an upward rigid shift, to the experimental 
energy value of the band gap [igI, of all the conduction band states. In general, for metals 
and conductor materials, the optical properties calculations demand a great number of 
eigenvalues and eigenvectors, therefore, a study of convergence, that considers the number 
of k-points used, is necessary. In the case of the insulating oxide Zr02 studied here, we have 
verified that no more than about 100 symmetrized k-points were sufficient. 

Figure 3 shows the imaginary part of the dielectric function (electronic contribution), for 
the three phases of Zr02, considering energies till 30 eV. For the cubic phase, all the diagonal 
components of 62 are identical (e^x = ^yy = ^zz)-, whereas the off-diagonal components are 
zero. The energy spectrum was adjusted by shifting it upwards, to the band gap experimental 



value for the monoclinic phase. Eg = 5.83 eV 16|. For tetragonal and monoclinic phases. 



the adjustment was taken by considering the same experimental value as for the monoclinic 



phase, once this later is the observed phase at ambient temperature. The function 62, for 
t-Zr02, has all off-diagonal components equal to zero, but symmetry assures that Sxx = 
Eyy 7^ Ezz- For the monoclinic phase, the components that are different from zero are e^xi 

^yyi ^zzi &nd Bxz- 

In order to avoid confusion, results for both full-relativistic (frel) and non-relativistic 
(nrel) calculations are shown in Fig. 3 only for the cubic phase. Relative changes in the 
values of the dielectric function, due to relativistic effects amount to 5%, at maximum, 
for all cases. For the tetragonal and monoclinic phases, the presented results derive from 
full-relativistic calculations. A comparison of the theoretical findings with the experimental 
result of Camagni et al. 2J], as derived from reflectivity data, shows better agreement for 
the monoclinic phase for which two more pronounced peaks occur in the theoretical curve, 
at 7.5 eV and ~ 11.5 eV. 

The real part of the complex dielectric function, Si, is presented in Fig. 4, for the three 
phases of Zr02. As for the case of 82 (shown in Fig. 3), we show results for Si as obtained from 
full-relativistic and non-relativistic calculations only for c-Zr02; for the other two phases, 
tetragonal and monoclinic, the results depicted in parts (b) and (c) of Fig.4 correspond to 
full-relativistic calculations. 

In order to compare the theoretical result with experiment, we show in Fig. 5 the calculated 
xx-component of the real part of the dielectric function, for m-Zr02, till energies of 
30 eV, together with the data reported by Camagni et al. 2S\ for two samples of m-Zr02 
stabilized with 12% and 24% of Yttria. Excellent accordance is seen between theory and 
experiment. 

In Fig. 6, we show the calculated reflectivity (R) for the three phases of Zr02, together 
with the experimental results, as extracted from French et al. [16!] and from Camagni et 
al. {24]. In the former work, the authors have considered three different samples of Zr02, 
in cubic, tetragonal, and monoclinic phases. These experimental results are referred as 
Exp. [2] in Fig.6. The data in the figure, referred as Exp.[l], correspond to m-Zr02 single 
crystals stabilized with 12% of Yttria 2^. Again, a good agreement between theoretical 
and experimental results is observed for the monoclinic phase. 

Figure 7 shows the behavior of the refraction index, as function of energy, for all three 
phases of Zr02, as derived from full-relativistic calculations. For the cubic phase, we also 
included the result obtained from a non-relativistic calculation, for comparison. In Table III 
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are depicted the values of the refraction index, at two different energies, for the tetragonal 
phase. They compare well with the data reported by French et al, as obtained from VUV 
spectroscopy measurements [l^. 



E. Phonon Contribution to the Dielectric Function 

In order to realistically describe the dielectric function and, in particular, to evaluate 
the static dielectric constant, £1(0), we have also to consider the phonon contribution to 
the low-frequency region of the frequency-dependent dielectric function. The enhancement, 
observed in the measured values of the Zr02 dielectric constant, is due to the role played by 



the optical phonons. We have used the Abinit code [22] to unveil this phonon contribution to 
the dielectric functions of the three phases. With this program, we calculated the equations 
of state, the phonon frequencies and eigendisplacements for the structures (not shown in 
this work). In this case, we adopted the Troullier-Martins pseudopotentials in the calcula- 



tions [25|, and the phonon dynamics were obtained by means of the adiabatic perturbation 
density functional theory. According to the available data, and previous theoretical results, 
for the phonon frequencies [14], our results are in complete agreement. 

From the obtained results, we have calculated the frequency dependence of the^real part 



of the dielectric constant (as described by Rignanese et al. 17| and Gonze et al. [26]) for the 
monoclinic, tetragonal, and cubic phases. 

In Fig. 8 we show the results for ei"'^^\w), for m-Zr02, along the x, y, and z axis, together 
with a fitting by means of the Lorentz-Drude model, i.e. 



e--' ^uj) = J2 u, ,2 , , ^.2, .21 ' (2) 



where n is the total number of oscillators, each one with the transversal optical (TO) phonon 
frequency utoi, ^pi is the plasmon frequency for the i^^ oscillator, and 7^ is its damping 
parameter, which is directly related to sample size. The Drude-Lorentz curves are a guide 
for the eye, and show that our calculated data follow this model. In the case of the x- 
component of ei"'*^\w), we included four oscillators, and the obtaining fitting parameters 
were Upi = 541.83 cm~^ and 71 = 0.05 cm^^ for the TO mode at 709 cm~^, ujp2 = 635.54 
cm~^ and 72 = 0.05 cm~^ for the TO mode at 485.1 cm~^, = 1052.34 cm~^ and 73 = 
0.05 cm~^ for the TO mode at 316 cm~^, and Up4 = 493.07 cm~^ and 74 = 0.05 cm~^ for 
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the TO mode at 220 cm~^. 

To describe the frequency dependence of the dielectric constant parallel to the y axis, given 
by the Lorentz-Drude model, we needed five oscillators. The obtained fitting parameters 
were Upi = 517.81 cm^^ and 71 = 0.05 cm~^ for the TO oscillator at 576 cm~^, ujp2 = 493.44 
cm^^ and 72 = 0.05 cm^^ for the TO oscillator at 473 cm~^, = 997.54 cm~^ and 73 = 
0.05 cm~^ for the TO oscillator at 413 cm~^, Up4 = 576.96 cm~^ and 74 = 0.05 cm~^ for the 
TO oscillator at 354 cm~^, and cUps = 465.39 cm"^ and 75 = 0.05 cm~^ for the TO oscillator 
at 247 cm~^ 

Finally, we have also to include five oscillators, in order to describe the frequency depen- 
dence of the dielectric constant parallel to the z axis. The obtained fitting parameters were 
Upi = 147.44 cm^^ and 71 = 1.22 cm~^ for the TO mode at 485.1 cm~^, Up2 = 353.93 cm~^ 
and 72 = 0.81 cm~^ for the TO mode at 399 cm~^, Ups = 1029.0 cm~^ and 73 = 2.18 cm~^ 
for the TO mode at 354 cm~^, Up^ = 513.52 cm~^ and 74 = 11.1 cm~^ for the TO mode at 
316 cm~^, and Up^ = 327.06 cm~^ and 75 = 0.0001 cm^^ for the TO mode at 220 cm^^. 

The calculated static dielectric tensor, due to lattice contribution, for the monoclinic 
phase of Zr02 is 



^ 24.12 1.77 ^ 
21.20 
1.77 17.02 



(3) 



Therefore, using £1(0) 
get 



_{electr) , 



(0) +£^""^(0), and using the e^i''^'^^^\Qi) from Figure 4(c) we 



(electr) , 



1 29.52 1.83 ^ 

26.20 . (4) 
\ 1.83 21.92 y 

The average value of the static dielectric constant, for the monoclinic phase, is then given 
by £1(0) ~ 25.9, which agrees fairly well with reported experimental values of ~ 25 [Ref.j2|]. 

It is worth mentioning that, we have also evaluated the dielectric constants for t- and 
m-Zr02, by using the same approach described above. The calculated phonon contribution 
for the tetragonal phase, added to the corresponding electronic part, lead to an average value 
of the static dielectric constant of £i(0) = 55.8, a value which is as high as that obtained 
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for the cubic phase, 44.4. Details of these calculations will be published elsewhere |27l |. 
These findings are another indication that the measured dielectric constant for the Zr02 is 
correctly described by considering the oxide in its monoclinic phase. 



IV. CONCLUSIONS 

We have performed an extended study, based on ab initio calculations, of the structural, 
electronic, and optical properties for the cubic, tetragonal, and monoclinic crystalline phases 
of Zr02, which takes into account scalar relativistic and spin-orbit contributions to the 
electronic structure, as well as the phonon contributions to the low-frequency regime of the 
dielectric function. For the band structure calculations, we used the DFT-GGA approach 
by means of the FLAPW method. Due to the spin-orbit interaction, we observe a spin-orbit 
(so) splitting energy of the top of the valence band (v) at the F-point, which for the cubic 
phase is Ag^(F) = 69 meV, while for the tetragonal phase is A^q(F) = 9 meV. So far, there 
are no experimental data reported for A^o in Zr02. 

The carrier effective masses are shown to be highly anisotropic, with the relativistic 
corrections playing an important role; this is of special interest when modelling advanced 
Zr02-based devices. In addition we report, with the relativistic corrections, an evaluation of 
a spin-orbit effective mass in the valence band, as well as the heavy- and light-hole (electron) 
effective masses at the valence (conduction) band, in particular for cubic Zr02. 

The results obtained for the real and imaginary parts of the dielectric function, for the re- 
flectivity, and for the refraction index, show good agreement with the available experimental 
results. In order to obtain the static dielectric constant of Zr02, we added to the electronic 
part, the optical phonons contribution, which leads to values of £i(0) ~ 29.5, 26.2, 21.9, along 
the XX, yy, and zz directions, respectively, for the monoclinic phase, in excellent accordance 
with experiment. 

Relativistic effects, taking also into account the spin-orbit interaction, are demonstrated 
to be important for a better evaluation of the effective mass values, and in the detailed 
structure of the frequency dependent complex dielectric function. We were able to unveil 
and estimate a more realistic total value of the static dielectric constant of Zr02, in the 
form £1(0) = ef'^'^\0) + e'f*^{0), yielding an average value of ~ 25.9, in good agreement with 
reported experimental values of ~ 25. 
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Figure 1 

FIG. 1: Band structures, along high symmetry axis of the BZ, and total density of states (TDOS) 
for (a) cubic, (b) tetragonal, and (c) monoclinic phases of Zr02, as obtained from non-relativistic 
calculations. The energy zero was taken at the valence band maximum, shown by a dashed hori- 
zontal line. The main character of the peaks in the TDOS is emphasized. 
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Figure 2 

FIG. 2: Band structures, along high symmetry axis of the BZ, and total density of states (TDOS) 
for (a) cubic, (b) tetragonal, and (c) monoclinic phases of Zr02, as obtained from fuU-relativistic 
calculations. The energy zero was taken at the valence band maximum, shown by a dashed hori- 
zontal line. The main character of the peaks in the TDOS is emphasized. 
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FIG. 3: Imaginary part of the dielectric function, versus energy, for the (a) cubic, (b) tetragonal, 
and (c) monoclinic phases of Zr02, as obtained from full-relativistic calculations. Only in (a), 
results as obtained from full-relativistic calculations (/re/) and non-relativistic calculations (nrel) 
are shown, for comparison. The theoretical curves were adjusted by performing a rigid shift upwards 
in energy to the experimental gap value for the monoclinic phase E^^P = 5.83 eV [Ref. [l^]. 
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Figure 4 

FIG. 4: Real part of the dielectric function, ei, versus energy, for the (a) cubic, (b) tetragonal, and 
(c) monoclinic phases of Zr02, as obtained from full-relativistic calculations. Only in (a), results as 
obtained from full-relativistic calculations (frel) and non-relativistic calculations {nrel) are shown, 
for comparison. The theoretical curves were adjusted by performing a rigid shift upwards in energy 
to the experimental gap value for the monoclinic phase Eg^^ = 5.83 eV [Ref. [if 
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FIG. 5: xx-component of the real part of the dielectric function, £ixx, for monoclinic Zr02, as ob- 
tained from full-relativistic FLAPW calculations, together with the experimental results as derived 
from reflectivity data, extracted from Ref.[24] for two yttria-stabilized samples. 
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FIG. 6: Reflectivity (in %) versus energy, for the (a) cubic, (b) tetragonal, and (c) monoclinic phases 
of Zr02, as obtained from fuU-relativistic calculations. Only in (a), results as obtained from full- 
relativistic calculations (/re/) and non-relativistic calculations {nrel) are shown, for comparison. 
Also shown are the experimental curves, as extracted from French et al. [16] (Exp.[l]) and from 
Camagni et al. B (Exp. [2]). 



21 




FIG. 7: Refraction index versus energy for the (a) cubic, (b) tetragonal, and (c) monoclinic phases 
of Zr02, as obtained from full-relativistic calculations. Only in (a), results as obtained from full- 
relativistic calculations {frel) and non-relativistic calculations (nrel) are shown, for comparison. 
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FIG. 8: Phonon contribution (full circles) to the low-frequency dependence of the real part of the 
dielectric function, as obtained from ab initio calculations, parallel to the x, y, and z 

axis, for monoclinic Zr02. The solid line corresponds to a fitting by means of the Lorentz-Drude 
expression, and was used as a guide to the eye. 
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TABLE I: Band gap energies, Eg (in eV), and symmetry of the valence-to-conduction band tran- 
sition, obtained from the non-relativistic {E^'^^) and full-relativistic {E^'^^^) calculations for the 
cubic, tetragonal, and monoclinic phases of Zr02. 

phase valence-to-conduction band transition symmetry E^*^^ E^'^^'- 

cubic X — 3.09 3.30 

r — ^ r 3.61 3.80 

X — > X 3.65 3.72 

tetragonal Z — >V 3.80 4.01 

X — ^ r 3.83 4.04 

r — ^ r 3.88 4.09 

monoclinic T — > X 3.44 3.58 

r — ^ r 3.82 3.98 

X — >X 3.50 3.64 
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TABLE II: Valence and conduction band effective masses (in units of the rest free electron mass, 
mo), at relevant symmetry points of the BZ, for the cubic (c), tetragonal (t), and monoclinic (m) 
phases of Zr02, as obtained from full-relativistic band structure calculations. In parenthesis, are 
shown the effective mass values as derived from non-relativistic calculations, m^, m*, m|g, 
m^^, m^*^, and m*^ stand for, respectively, hole, electron, heavy-electron, light-electron, heavy-hole, 
light-hole, and split-off hole effective masses. 



phase direction valence conduction 

c r — > L m*^ = 0.26;mf^ = 0.23(0.23);m*^ = 0.28(0.27) ml^ = 0.28(0.26);^;*^ = 0.27(0.26) 

r — > X m*o = 0.77;m^f^ = 0.36(0.24);m^^ = 3.84(4.23) ml^ = 2.09(2.06);m;*g = 0.51(0.48) 

X — ^ r ml = 0.28(0.27) m* = 1.77(1.97) 

X — >W ml = 3.38(3.29) m* = 1.17(1.33) 

t M^r = 0.47(0.47) ...(...) 

r — > Z ml = 5.31(6.68) m* = 2.32(2.32) 

Z — >R m* = 0.83 (0.88) ...(...) 
m X — ...(...) m* = 1.16(1.12) 

X — > S ml = 1.09(1.06) ml = 2.32(2.32) 

r — >Y ml = 2.55(3.08) m* = 1.39(0.75) 

r — > Z ml = 1.49(1.42) m* = 4.08(4.18) 



TABLE III: Values for the refraction index, at certain energies, as obtained from non-relativistic 
(n"*"*^') and full-relativistic (n-^*"^') calculations, for tetragonal Zr02. For comparison, the experi- 



mental results as observed from VUV spectroscopy measurements, and extracted from Ref. 
are shown in the last column {n^xp)- n^x and Uzz stand, respectively, for perpendicular and parallel 
refraction index to the axis along c, one of the lattice constants in the tetragonal lattice. 
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